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Abstract 

We calculate the one-loop supersymmetric electroweak-like and QCD-like 
corrections to the top-quark pair-production cross section at the Tevatron, in- 
cluding the important effects of non-degenerate top-squarks and left-right top- 
squark mixing. The largest electroweak-like effects yield a negative shift in the 
cross section and are enhanced right below the threshold for top-quark decay 
into top-squark and higgsino-like neutralino, and can be as large as —35%. 
The largest QCD-like effects are positive and are enhanced for light top-squark 
masses, and can be as large as 20%. Such shifts greatly exceed the present 
theoretical uncertainty in the Standard Model prediction, and therefore may 
be experimentally observable. We also explore the one-loop shifts in scenar- 
ios containing light top-squarks and higgsino-like neutralinos that have been 
recently proposed to explain various apparent experimental anomalies. 
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1 Introduction 



The existence of the top quark has now been firmly established by the CDF and 
DO Collaborations at Fermilab ||T|. The analysis of over lOOpb -1 of data by each 
experiment has allowed a determination of the top-quark mass and pair-production 
cross section || 

mf DF = 176 ± 9GeV, ^ = 7.61^ P b ; (1) 
m™ = 170±18GeV, af { ° =5.2±1.8pb; (2) 

to better than ~ 5% and ~ 20% respectively. The experimental error flags are 
expected to be further reduced when the data set is fully analyzed. A significant 
increase in sensitivity will become available once the Main Injector upgrade becomes 
operational in 1999, with an expected reduction in the top-quark mass uncertainty 
down to 3.5 (2.0) GeV and in the production cross section down to 11% (6%) with 
an integrated luminosity of C = 1 (10) fb _1 ||. On the other hand, the theoretical 
prediction for the cross section in the Standard Model has become rather precise 
0, [|, ||, and is presently known to better than 10% for fixed values of m t , e.g. || 

^theory ^ = g ^+0.09 pb ? ^theory ^ = g ^+0.07 pb ( g } 

The agreement between theoretical expectations and experimental observations is not 
tight enough to preclude moderate shifts (~ 2 pb) to the production cross section from 
new physics phenomena. However, greater than 50% shifts in the Standard Model 
prediction due to new physics may conflict with experimental observations. It is then 
opportune to investigate how large a shift scenarios for new physics may yield, most 
interestingly in the case of low-energy supersymmetry. 

In this paper we consider the one-loop corrections to the pair-production cross 
section for top quarks at the Tevatron, including both supersymmetric electroweak- 
like (loops of charginos and bottom-squarks or neutralinos and top-squarks) and 
QCD-like (loops of gluinos and squarks) contributions. Our calculation extends and 
corrects those in the literature J7], [|, identifies the largest contributions and deter- 
mines under what conditions they may be enhanced, and treats all important and 
realistic effects (such as top-squark mass non- degeneracy and mixing) in a unified way. 
We also consider specific scenarios with especially light supersymmetric particles that 
have been proposed in connection with the Rb observable |§, an as-yet unexcluded 
light higgsino window Jl0|, and a possible supersymmetric explanation Jy], [12] for 



the ecyy event observed by CDF [[y|. We stress that light supersymmetric particles 
that affect the top-quark cross section are also likely to affect its decay width via new 
supersymmetric channels, and the resulting decrease in B(t — > bW) may affect the 
yield of bWbW events in a much more significant way than a shift in the underlying 
production cross section. 

This paper is organized as follows. In Sec. we present compact expressions for 
the results of our analytical calculations for the one-loop supersymmetric electroweak- 
like and QCD-like corrections and contrast them with existing calculations. In Sec. [3] 
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we present a qualitative and quantitative study of the expected corrections and iden- 
tify regions in parameter space that entail the largest possible effects. These effects 
may easily exceed the theoretical uncertainty in the Standard Model prediction for 
the top-quark cross section, and thus may be disentangled at future Tevatron runs. In 
Sec. £| we adapt our general results to study specific models that have been proposed 
in the literature and that contain especially light supersymmetric particles. Sec. |] 
summarizes our conclusions and the Appendix contains some further expressions. 



2 Analytical results 

Top-quark pair-production at a hadron collider proceeds at tree-level via two under- 
lying s-channel parton processes: qq annihilation (see Fig. 0(a)) and gluon fusion. At 
the Tevatron the first process dominates (90%) M, and therefore we neglect the latter 
in what follows. One-loop corrections to qq — > g — > tt in the context of low-energy 
supersymmetry fall into two categories: electroweak-like |7| and QCD-like ||. The 
significant electroweak-like vertex and external leg corrections modify the outgoing 
part of the Feynman diagram via new contributions to the gti vertex and to the 
normalization of the top-quark wavefunction that involve loops of top-squarks and 
neutralinos or bottom-squarks and charginos (see Fig. The QCD-like corrections 
modify both the incoming and outgoing portions of the diagram in an analogous 
fashion, but involve loops of gluinos and first-generation-squarks or gluinos and top- 
squarks (see Fig. |]). There are also (smaller) QCD-like box diagram contributions 
involving gluinos and squarks (see Fig. |3](i)). 

The tree-level qq — > tt parton cross section is given by 

*=^/M(? + 2rn?) (4) 



where /3 t = y 1 — 4mf/s, and Vd is the usual parton- level center-of-mass energy. The 
supersymmetric electroweak-like one-loop corrections involving the Higgsino compo- 
nents of the neutralinos and charginos are enhanced by a large top-quark Yukawa 
coupling as seen below 

At (^j H 2 t R 3 At h L R ± t R , X t t L Hl t R , \ t TlHI t R , (5) 

where the carets represent superfields and the terms on the right-hand side include 
the Yukawa coupling components of interest. These corrections involve loops of top- 
squarks and neutralinos or bottom-squarks and charginos and are given by 



v — .- . . — / v \ 

EW 




^+2m < 2 )(F™ + F 1 c ) + 2(F 5 n + F 5 c )m t s , (6) 



where X t = gm t / (V%M W sin/5) is the top-quark Yukawa coupling, and tan/5 = v 2 /vi 
where vi,v 2 are the Higgs vacuum expectation values that arise in the MSSM. Also, 
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F"g C are form factors that encode the loop functions and depend on the various mass 
parameters. The top-squark-neutralino form factors are given by: 



(X p tj,tj) 



F? 



J2\N jA N; 4 £ [c 24 + m 2 t (c n + c 21 ) + \B X + m 2 t B[ 

j=l L J=l 

-Re (N j4 N j4 ) £(-l) J+1 M^t) m t m x o (c + c n + B' Q )W> ij >^ 
j=i 

E{^iV; 4 £ [-\m t { Cll + c 21 
j=l L J=l 



(7) 



(X°,tj,tj) 



+Re (iV j4 JV i4 ) E (-1) J+1 1 sin(20 t ) m x o (c + c n 



j=i 



Here iVj 4 represents the higgsino admixture of x°j {H® m ^ ne notation of Eq. (|5|) and 
Ref. [0), and the J = 1, 2 sum runs over the two top-squark mass eigenstates (m^ 2 ), 
which are obtained from the gauge eigenstates via: t\ = cos9 t tt + sin 6 t tR and 
ti = — sin 6 t ti + cos #t £r. The various B and c functions in the above expressions are 
the well documented Passarino-Veltman functions Jl5[ (adapted to our metric where 
); the B functions depend on (m t , m x o, m^) whereas the c functions depend 



P? 



on (— P3,P3 +P4, iTT'x ., m lji m tj) [ as reminded by the superscripts in Eqs. where 
P3 and pi are the momenta of the outgoing top-quark and anti-top-quark respectively 
(see Fig. |lj(b) ) . The bottom-squark-chargino form factors are in turn given by 




( cos 2 


or 




\ sin 2 


ov- 




/ cos 2 


op 




\ sin 2 


»r 





c 2 a + m 2 (cn + c 2 i) + \B X + m 2 B' l 



(Xf,bj,h 



jm t (cn + c 2 i, 



(Xj,bj,bj) 



9) 
(10) 



where Vj 2 represents the higgsino admixture of xtt an d f° r completeness we have 
allowed a non- vanishing bottom-squark mixing angle such that b\ = cos 6 b 6 L +sin 6 b bn 
and 6 2 = — sin b\ + cos 9b bR. The J = 1,2 sum runs over the two bottom-squark 
mass eigenstates (rag ). Note that if ^ = 0, only the left-handed bottom-squark 
is involved in the loops. In Eqs. (p|JT0|) the B functions depend on (mt,m ±,mg J ) 
whereas the c functions depend on (— P3,p% + p 4 ,m ±,mg , rag ). 

The functions B\ and c 24 in Eqs. (0,fD contain infinities. We used the modified 
minimal subtraction scheme and introduced counterterms to eliminate them. Never- 
theless, it turns out that the infinities (and the renormalization-scale fi dependence) 
of these two functions cancel each other out in both equations (independently) even 
without introducing counterterms. 

Since the Passarino-Veltman functions can be notoriously difficult to evaluate 
numerically for certain values of the parameters, we have employed different methods 
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(all of which agree), in particular using the software package ff JT6|. Our results 
for the case of no top- or bottom-squark mixing (9 t = Q b = 0) disagree with those 
originally published in Ref. However, these authors have since revised their cal- 



culation and their corrected expressions [T7| are now in agreement with our results 
above. Our results for the realistic case of squark mixing (9 t , 0& ^ 0) are new.Q 

In Eqs. (|6||T||8l,P,|TCil) we have not explicitly exhibited the additional contribu- 
tions to the Fx§ form factors that arise from the gaugino admixtures of the neutralinos 
and charginos. These are proportional to the electroweak gauge couplings and are 
therefore not enhanced by a large top-quark Yukawa coupling. These contributions 
have been neglected in our present analysis. For completeness, in the Appendix we 
present analytic expressions for the contributions from all (higgsino and gaugino) ad- 
mixtures to the vertex and external leg corrections. We have also not exhibited the 
electroweak-like corrections to the incoming part of the diagram, that involve loops 
of first-generation squarks and neutralinos or charginos. These corrections are rather 
small as the relation m q pa m q + m x , that the numerical analysis below shows is 
required for enhancement, is not satisfied in this case. Furthermore, for higgsino-like 
neutralinos these contributions are proportional to m 2 (c.f. Eq. (|6|)) and are therefore 
negligible. 

The Q CD-like corrections involve self-energy, vertex, and box diagrams. In the 
approximation of neglecting the box diagrams,^] one can readily obtain the shifts in the 
total parton-level cross sections and cast them in the same format as the electroweak- 
like corrections given above. In this form comparison between the two classes of 
corrections becomes transparent. The QCD-like one-loop correction becomes 

A(?Q cd = ^SprA (g) [§(? + 2m?)(*? + F?) + 2F 5 W] • (11) 
The form factor F{, which modifies the outgoing part of the diagram is given by 

F{ = E{(-|)[c 24 + m?(c 11 + c 21 )] ( ^ ) + (|)[i?i + 2mM 
j=i 

+(-l) J sm(2e t )m t m~ g [(l)(c + c u )^ - 
+ (f)[ _ I + 2c 24 + s(c 22 - c 23 ) - m 2 g c - m 2 (c + 2cn + c 2 i) 

-(-l) J 2 S m(26 t )m t m~ g (c + Cn)f J§§) } , (12) 

where 6 t is the top-squark mixing angle, the B functions depend on (m t , m g , m^), 
and the superscripts (abb) indicate that the corresponding c functions depend on 
(~P3,P3 + P4,m a ,m b ,m b ). 

The analogous form factor that modifies the incoming part of the diagram (Ff) 
is obtained from F[ by setting ij — > qj, m t — > m q , 9 t — > 9 q , and by replacing ps,p4 

1 While the present paper was being written up Ref. appeared, which contains a calculation 
of the 0t ^ case, that agrees with our result above in the case of real Nj4 values. 

2 Recent calculations of supersymmetric QCD-likc corrections to light-quark scattering at the 
Tevatron indicate that box diagram contributions are indeed small near the threshold region [|19| . 
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by Pi,P2, the incoming quark and anti-quark momenta. Since m q « 0, the expression 
for Fi simplifies considerably. Note the close resemblance between the first two lines 
in Fl and the expression for F™ in Eq. ([7|) above, as both form factors originate from 
analogous diagrams (jq «-> g); the differences in the coefficients stem from the color 
factors that appear in the QCD-like diagram. In contrast, the last two lines in F{ 
correspond to a new diagram (not present in the electroweak-like case) where the 
gluino couples directly to the gluon and which carries a large color factor (Fig. |||(g)). 
We also have 



F i = E{(|)K( Cll + c 21 )-(-l) J sin(2^K(c + Cll )] 
j=i 



(gtjtj) 



+ (|)K(c n + c 21 ) + (-l) J sin(2^)m^c 11 ] ( ^ ) } , (13) 

where we again note the resemblance (up to the color factor coefficient) between the 
first term in F§ and its counterpart F§ in Eq. (H). The second term in F§ corresponds 
to the diagram not present in the electroweak-like case. The analogous form factor 
Ft- (that arises from the incoming part of the diagram) is obtained from F§ by setting 
tj — > qj, m t — > m g , 9 t — ► 6 q , by introducing an overall minus sign and by replacing 
P3)P4 by pi,P2, the incoming quark and anti-quark momenta. This contribution is 
negligible since m q ~ and because in the MSSM 9 q oc m q « 0. Our results for the 
super symmetric QCD-like corrections agree with those presented earlier in Ref. ||, 
and its erratum [fJU] . 



For general neutralino composition there exist box diagrams that mix electroweak- 
like and QCD-like corrections,^ with gluinos, squarks, top-squarks, and neutralinos 
in the loop (see Fig. |4j). For the higgsino-like neutralino case that we consider, these 
diagrams are proportional to the light-quark Yukawa coupling and therefore negligi- 
ble. 

The actual observable cross section is obtained by integrating a + Act over the 
parton distribution functions, i.e., 

cr + Act = / dr / m(xi)m(x 2 ) + u(xi)sea(x 2 ) + sea(xi)u(x 2 ) 

J TO J T X\ 

+d(xi)d(x2) + d(xi)sea(x2) + sea(xi)d(x2) 

+6sea(xi)sea(x 2 ) (a + Aa)(s) , (14) 

where r = imf/s, x 2 = t/xi, s — rs, and the parton distribution functions (u, d, sea) 
are taken from Ref. [21] setting the scale Q = m t . 



3 We thank C. Kao for point this out to us. 
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3 Numerical results 



Inspecting the above formulas, one can immediately get an idea of the typical size of 
the one-loop corrections, as they are proportional to 

— K w ' — K UJ ■ (15) 

For the favored values of m t , both these factors are ~ 1%. We thus see that unless 
there are enhancements within the loop factors contributing to these diagrams, both 
contributions are comparable in size, and more importantly, much too small to be 
disentangled from the Standard Model contribution, or to be observed experimentally 
at the Tevatron or any of its planned or proposed upgrades. 

Dynamical enhancements are however possible for restricted ranges of the mass 
parameters. The top-quark self-energy diagrams show one such enhancement through 
the B' ,B[ functions when m t ~ m a + mj, where m a ^ are the masses of the two 
particles entering the two-point function. The enhancement occurs right below the 
threshold for t — > a + b decay. For the top quark there are three possibilities for m a} b'- 
(m^ 2 , m x o 2 3 4 ) ; (jt^ 2 , m x ± 2 ) ; and (m^ 2 , uig) . Given the present experimental lower 

limits on the squark (excluding t) and gluino masses (i.e., mg,m g > 175 GeV; nig rs 
rrig > 230 GeV p2||), only the first possibility may be realized. That is, enhancements 



may occur for m, m mx „ +m v o . Moreover, these enhancements will be maximized 

'liS Xl, 2,3,4 

when the neutralinos have a high higgsino content. (We do not see such enhancements 
for the incoming part of the diagram because m q ^ 0.) 

In this region of parameter space (and for real values of Nj±) one obtains the 
following approximate expression 

Acr EW / A 4 x 2 4 



V / j=l J=l 

^ (^) p2[m 2 t B[-(-iy^sixi(2e t )m t m x B' 



(16) 



where the second expression follows when one of the neutralinos (x) carries the full 
higgsino admixture. Also, the choice tan/3 = 1 maximizes A t for a fixed value of m t . 
In this 'best case scenario', one can plot Aa/a (i.e., after integration over parton 
distribution functions) versus the neutralino mass (m x ), and study the dependence 
on the top-squark masses and mixing angle. Note that the mixing-angle term does 
not contribute if the top-squark masses are degenerate. 

In the case of degenerate masses, taken at = = 50 (75) GeV, the 
resulting relative shift (%), as a function of m x , is shown by the dotted curve on 
the upper-left-hand panel in Figs. ||(|D.f3 The large dips at m x w m t — m^ 12 ~ 

4 The numerical results in Figs. |B| and ^ have been obtained for the 'best case scenario' outlined 
below Eq. (|l6|) and include the complete electroweak-like vertex and external leg corrections. 
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125 (100) GeV derive from the B[ term in Eq. (|T6| ) and have been regularized by 
setting m\ — > m\ — im t T t , where T t (a few GeV) is the top-quark decay width|] The 
depth of the dips depends on the value of m x /m t at which it occurs, and is maximized 
for m x /m t ~ 0.73, which in this case implies m x pa 128 GeV, m 4 ~ pa 47 GeV (i.e., as 
in Fig. |5] (dotted curve, upper- left-hand panel)). 

In a more realistic scenario, light top-squarks cannot be degenerate in mass 
because limits on additional contributions to the p parameter restrict the splitting 
between the and bi ~ bi, and from direct experimenal searches one estimates that 
mr > 200 GeV. In the case of no top-squark mixing (6 t = 0, |), the splitting of tip 
leads to a double-dip structure if both top-squark masses are such that m t pa +m x 
can be satisfied, as in Fig. [| (solid curve, upper-left-hand panel), where we have taken 
= 50 GeV and m 4 ~ 2 = 100 GeV. Note that the dips do not have the same depth, as 
discussed above. In Fig. |6] there is a single dip because we have taken = 75 GeV 
and m h = 250 GeV. 

Yet more realistic is the case of top-squark mixing, which is naturally present 
in supergravity theories. This mixing tends to "screen" the contributions of top- 
squarks to the Z-pole observables p3| , and therefore allows a larger ti — &l mass 



splitting. In Figs. [| and |6] we present the results for four choices of the mixing angle 
9 t = 0,0.10,0.25, and maximal mixing (|). The plots also apply for d% — > | — 9 t , 
which leaves sin^t) unchanged. In this case the B' function in Eq. (|TIJ) plays an 
important role, as it exhibits a similar dip behavior as B[ does, although with a 
different sign at each dip. This effect makes one dip deeper (corresponding to ti) 
whereas the other one shallower (corresponding to t 2 ) ? as evidenced in Fig. |5|. The 
effect can be very significant, completely eliminating one of the dips at maximum 
mixing angle. In Fig. |5] only the dip that gets deeper exists. 

As explained above, the QCD-like corrections are not expected to exhibit the 
dip structure that the electroweak-like corrections possess, because the relation m t ~ 
+ m- t cannot be satisfied by the experimentally allowed gluino masses. In Fig. |7] 



ni 



we show the QCD-like correction versus the universal squark masses (mi 2 = m^) for 
fixed values of the gluino mass. The curve for m g = 150 GeV behaves differently from 
the others because for sufficiently light squark masses, the relation m t pa m~ g + m g 
will be satisfied, i.e., a dip occurs. This region of parameter space is disfavored 



experimentally; it is shown here to make contact with Ref. pOJ, with which we agree 
qualitatively. From Fig. [7] we see that shifts as large as ~ 20% may occur. We 
also observe that the corrections go to zero when the supersymmetric particle masses 
get large, i.e., the expected decoupling effect. This figure also contains a curve (the 
dashed line) where all sparticle masses are taken to be the same (mg = m% x 2 = m^), 
which is seen to intercept the other curves at the appropriate places, and to decouple 
rather quickly. 

The significant QCD-like shifts observed in Fig. [7], especially for m~ g = 200 GeV, 



5 Throughout our numerical calculations we have used mt = 175 GeV and m t r t = 289 GeV 2 . For 
a s we took the world-average value of 0.118 at the scale /i = Mz in the modified minimal subtraction 
scheme. 
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may be understood in terms of the large color factor in the vertex correction diagrams 
that couple the gluon to two gluinos, both in the incoming (Fig. §(c)) and outgoing 
(Fig. 0(g)) parts of the diagram. A relevant role may also be played by a dynamical 
enhancement of the vertex diagrams that occurs for = 2m g , Irriq. As \/d is inte- 
grated over the parton distribution functions, starting at = 2mt and peaking at 
ps 400 GeV, the typical cusp is smeared out but it appears visible as the peak in 
the dashed curve in Fig. [F] (obtained for rriq = rrig). 

In order to explore the effects of lighter top-squark masses, we concentrate on 
the rriq = m g caseQ (the dashed line in Fig. [7]) and in Fig. |8| we plot the QCD-like 
corrections for representative choices of (m^, mz). For reference, the all-equal-masses 
case is shown as a dashed line (as in Fig. |7|). We can see that light top-squark masses 
enhance the corrections, especially in the (200 - 250) GeV. This 

enhancement occurs (although to a lesser extent) even if only one of the top-squarks 
is light. We have also explored the effect of top-squark mixing, which is no n- vanishing 
only if mi ^ m i 2 - We find this effect to be rather small in the case of the QCD-like 
corrections amounting, for example, to a decrease in the peak value of the (50,250) 
curve in Fig. |8| by 15% for maximal mixing. 



4 Expectations in specific models 

The results presented in the previous section should represent the largest one-loop 
super symmetric shifts (due to vertex and external leg corrections) to be expected in 
the top-quark cross section. However, in specific regions of MSSM parameter space 
or specific supergravity models, the shifts are likely to be much smaller than the 
largest possible ones, as the conditions for enhancement may not be satisfied: large 
electroweak-like shifts require a higgsino-like neutralino, light top-squarks, tan/3 ~ 
1, and a specific relation between their masses (i.e., m x + mi ~ mt); large QCD- 
like shifts require m 9 -, m g < 250 GeV and light top-squarks. Interestingly enough, 
various scenarios (i.e., selected regions of MSSM parameter space) recently proposed 
to possibly explain some experimental measurements that appear to deviate from 
Standard Model expectations, fall precisely in the class of models that may lead to 
one-loop enhancements of the top-quark cross section. 

The discrepancy between the LEP measured value of Rb and its prediction in 
the Standard Model may be alleviated by supersymmetric loop corrections to the Zbb 
vertex that involve charginos and top-squarks @. Moreover, tan/3 should be close 
to 1, the charginos should be higgsino-like (and correspondingly the neutralinos too), 
and the top-squarks should be right-handed. Both should be as light as LEP 1.5 
searches allow ||24j| . Since in this scenario there are no restrictions on the squark 



or gluino masses, let us concentrate on the electroweak-like corrections, that will be 
significantly enhanced (negatively) in this case if the relation m x + mi ps m t happens 
to be satisfied. This may indeed occur for mi 2 (but not for the desired values of m~ tl 

6 The rriq w rrig relation occurs naturally in supergravity theories. 
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and m x , both below Mw). 

In analogy with the so-called light-gluino window, it has been remarked that 
there is a light-higgsino window |i~0|| , where the Higgsino mixing parameter and the 
SU(2) gaugino mass are very small (//, M2 ~ 0) and tan/3 ~ 1. In this scenario there 
are three neutralinos with significant higgsino admixtures, two with masses close to 
Mz and a very light one. Moreover, a light top-squark is also desired to enhance R s b usy . 
Therefore, if m t « m x + m^ m Mz + mi or equivalently ~ m t — Mz ~ 85 GeV, then 
the top-quark cross section will be shifted to lower values by a significant amount. 



The last scenario we address has been advanced in Refs. |IT], |12f as a possible 
explanation for the one much-publicized event at CDF, consisting of ee'yy plus miss- 
ing energy. This event has been ascribed to selectron pair-production, with decay 
into electron and second-to-lightest neutralino an d further radiative decay of 

the neutralino into the LSP (x° ~~ * Xi + 7)- The missing energy is carried away by 
the pair of lightest neutralinos produced. The dominance of the (one-loop) radia- 
tive decay X2 ~> Xi + 7 over the more traditional ones (x° ~~ > Xiff) provides the 
most important constraint on the parameter space, requiring a higgsino-like Xi an d 
a photino-like x°- This is achieved by setting the SU(2) and U(l) gaugino masses 
equal at the electroweak scale (Mi = M 2 and tan/3 « 1), precluding the gaugino 
mass unification in GUTs. Furthermore, the kinematics of the event appear to re- 
quire: m x o « (30 — 55) GeV and m x o « m x o + 30 GeV. To make contact with the 
supersymmetric enhancement of Rb, it is also assumed that the top-squark is light 
(m^ « (45 — 60) GeV). In this case, however, the top-quark would have enhanced de- 
cays to t + Xi, thus diluting the observed top-quark sample. To undo this effect, it has 
been further proposed |T2| that squarks and gluinos should be as light as experimen- 
tally allowed (m g « (210-235) GeV, m- q « (220-250) GeV) such that g -> tt decays 
add to the top-quark sample significantly. (This mechanism was originally proposed 



in Ref. p5| .) Given these values of the squark and gluino masses, from Fig. [8] one can 
see that we generally expect m +15% supersymmetric QCD-like corrections to the 
top-quark cross section. (There may also be significant electroweak-like corrections if 
777£ 2 « m t — m x o « (105 — 145) GeV.) 

Going beyond the specific models discussed above, the light top-squarks that 
enhance both electroweak-like and QCD-like corrections will likely decrease the canon- 
ical top-quark branching ratio into bW because of the availability of the supersym- 
metric decay channels t — > t + x- Indeed, the top-quark Yukawa coupling in Eq. (|5|) 
entails an enhanced coupling between the top-quark, a higgsino-like neutralino, and 
a right-handed top-squark. This situation is favored by the enhanced corrections dis- 
cussed above, and therefore enhance the exotic decays of the top quark. At present 
there are no real experimental limits on B(t —> bW), only on B(t —>■ bW)/B(t — > qW) 



26]. The only limits on B(t -> bW) have been obtained by correlating the top-quark 
mass and cross section measurements with the Standard Model cross section, imply- 
ing B(t — > other) < 25% |27| . If the electroweak-like (negative) correction occurs 
when the exotic channel (t — > tx) is kinematically allowed, the yield of bWbW events 
will be decreased by the two effects, implying an effective B(t — > bW) ratio as small as 
(0.65) (0.5) ~ 0.3, assuming a —35% shift in the cross section and B(t — > tx) = 1/2. 
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However, we note that the enhancements to the electroweak-like corrections occur 
right below the threshold for top-quark decay into top-squark and neutralino (see 
Figs. |5|,|6|), and therefore the exotic decay channel is effectively closed. 

Recent complementary studies of top-quark properties in supersymmetric the- 
ories include supersymmetric one-loop corrections to the t t% exotic decay channel 
[ p8| , supersymmetric three-body decays of the top quark P5| , and one-loop super- 
symmetric corrections to the top-quark width |K| . 

5 Conclusions 

The experimental study of the top quark has just begun in earnest, with its mass and 
cross section having been measured to some precision. The Main Injector upgrade of 
the Tevatron should essentially provide a top-quark factory, where observations will be 
confronted with theoretical expectations for cross sections and branching ratios. We 
have shown that these precise measurements may indeed point to deviations from the 
Standard Model, as may be expected in supersymmetric theories. These deviations 
may be quite sizeable and therefore easy to detect, especially in scenarios with rather 
light sparticles that have been proposed to explain various apparent experimental 
anomalies. In the longer term, the presence of top quarks at the LHC will constitute 
one of largest backgrounds in new physics searches. Therefore, it will be essential 
to have a very good understanding of top-quark physics beforehand, so that these 
backgrounds may be subtracted off appropriately. 
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Appendix 

In this Appendix, we present analytical expressions for the one-loop supersymmet- 
ric electroweak-like vertex and external leg corrections to top-quark production by 
including all components of the neutralino (viz., the higgsino, the photino and the 
zino) and of the chargino (viz., the higgsino and the wino), as well as top-squark 
and bottom-squark mixing. Though the contributions are enhanced for higgsino-like 
neutralinos and charginos, in a detailed numerical calculation one would also need 
to include the contributions from gaugino-like neutralinos and charginos. We again 
ignore corrections arising from box diagrams. 
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The invariant amplitude for top-quark production via qq annihilation can be 
written as 

M = M + 6M (17) 

where M Q is the tree level amplitude and SM is the first-order electroweak-like cor- 
rection. M and SM are given by 



iM 
i5M 



v(p2)(-ig s T A Y)u(p 



u{pz){-ig s T A r)v{p4 



v{p2){-ig s T A Y)u{pi) Z f" M u (p 3 ) A^v (p 4 



(18) 
(19) 



where Pi,P2, Pz an d Pa are the incoming quark and anti-quark and outgoing top-quark 
and anti-top-quark momenta respectively. A M can be written as 

A" = -ig s T A (F^ + F 27 ^ 7 5 + F z k» + F 4 F 75 + F B ik v a^ + F e ik u a^ l5 ) (20) 

where k^ = p± + p% = p% +p% and a^ u = f [ 7 M , 7"]. The form factors i*i,..,5 encode the 
loop functions and depend on the various masses in the theory. The supersymmet- 
ric electroweak-like vertex and external leg corrections to the tree-level parton cross 
section are give by 



A? EW = ^ A [l(s + 2mf)(F? + F^) + 2(F 5 " + F 5 c )m t s 



(21) 



We see that the corrections depend only on F\ and F B . Below we present the complete 
form factors that appear in these equations. The tree-level parton cross section is 
given in Eq. ([|). 

The Feynman rules used to obtain the following form factors are given in 
Ref. [III]. For the electroweak-like corrections involving neutralinos and top-squarks 
we get 



cos 9 t + B n sin 9 t \ 2 + \C n cos 6* — A* sin 9 t 



x 



C24 + rn 2 (cn + c 2 i) + \Bi + m 2 t B[ 
-Re [(A n cos 9 t + B n sin 9 t )(C* cos 9 t - A n sin 9 t )\ 2m t m x o (c + c n + ff Q )W&&) 



-(ti <r-> t 2 , cos9 t <-» — sin# t , sin# t <-> cos 9 



t > 







j- 2 53 j (l^ncos^ + 5 n sin0 t | 2 + |C n cos#t - A* sin0 t | 



?mt(c u + c 2 i 



(X?,tl 



+Re [(A n cos 6> t + _B n sin 9 t )(C* cos 6*4 — v4 n sin m x o (cq + cn 



,(x?,ti,ti) 



-(t\ <-> cos^ <-> — sin# t , sin# 4 <-> COS6 1 
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where (in the notation of Ref. [[14]]) 



A n 
B n 



gm t 



2 s/2M w sin (3 



n* 4 = - 1 -\n* 4 



V2 



I 

72 



2 2sin 2 w 
- e N'* — 

3 3 3 cos 6w 
2 sin 2 



3 J 3 cos tV 



^; 2 



i 



2cos# 



-0^ 



J2 



(24) 
(25) 
(26) 



Recall also that t\ = cos 9 t ti + sin 9 t tu and £2 = — sin 9tti + cos 0$ £r. The various B 
and c functions in the above expressions are the well documented Passarino-Veltman 



functions |15] (adapted to our metric where pj 



ra- 



the B functions depend on 



(m t , m x q,mi 13 ) whereas the c functions depend on (— p s , p 3 + p 4 , m x q, 2 ,rn^ 1 2 ) [as 
reminded by the superscripts in Eqs. (22,23)]. We note that in the expressions for 
F" 5 given in the main text in Eqs. (0,H), which are specific to the case N'^ = iVj 2 = 0, 
we have extracted an overall factor of (Ai/4-7r) 2 . 

The form factors for the electroweak-like corrections due to loops involving 
charginos and bottom-squarks are given by 



cos^l 2 + \B c cos9 b + C c sin6> b | 2 ) 



X 



c 24 + rnf(c u + c 2 i) + \B X + m 2 t B' l 
-Re [A c cos 9 b (B* cos 9 b + C* sin 9b)} 2m t m ± (co + c\\ 



+ (6i 



bo, cos 9 h 



sin 9h, sin 9i, 



cos 9, 



n 



— Y\(\A C cos 9 b \ 2 + \B c cos9 b + C c sin 9, 



b) , 

2\ 



(27) 



-^m t (cn +c 2 i 



(xf,6i,6i) 



+Re [A c cos 9 b (B* cos 6 + C* sin b )] m ± (c + c u ) 



+(6i 



6 2 , cos6> fe 



sin^f,, sin^ 



COS ^ft 



(xfMM) 



where (in the notation of Ref. [|1~4 

A c - 

B c - 

a. -- 



2V2M w sin/3 j2 ' 
2^M w cos(3 j2 



W* 2 



-X b U j2 . 



(28) 

(29) 

(30) 
(31) 



Also 2 9b is the bottom-squark mixing angle defined such that b\ = cos 9 b &l + sin 9 b bR 
and b 2 = — sin 9 b b L + cos 9 b Sr. In this case the B functions depend on (m t , m x ±, J 
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whereas the c functions depend on (—^3,^3 + Pa, m x t i m b 12 ) m bi 2)' Note a ^ so that m 
the expressions for F° 5 given in the main text in Eqs. (PJlOD, which are specific to 
the case Uji = Uj2 = 0, we have extracted an overall factor of (A t /47r) 2 . 
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Figure 1: Tree-level diagram (a) describing qq — > tt production. The blob diagram 
(b) indicates the choice of external momenta used throughout our calculations. 
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Figure 2: Feynman diagrams describing the one- loop supersymmetric electroweak-like 
corrections to qq — > ti, including external leg corrections (a,b) and vertex corrections 
(c) from top-squark and neutralino loops. An analogous set of diagrams exists where 
the top-squarks are replaced by bottom-squarks and the neutralinos by charginos. 
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Figure 3: Feynman diagrams describing the one- loop supersymmetric QCD-like cor- 
rections to qq — > tt, including external leg corrections (a,b,e,f), vertex corrections 
(c,d,g,h), and a representative box diagram (i). 
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Figure 4: Representative box diagram (one of four possible ones) describing the mixed 
one-loop supersymmetric QCD-like and electroweak-like corrections. 
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Figure 5: The relative (%) one-loop supersymmetric electroweak-like correction to the 
top-quark pair-production cross section at the Tevatron as a function of the (higgsino- 
like) neutralino mass, for m t = 175 GeV, tan/5 = 1, = 50 GeV, m^ 2 = 100 GeV, 
and various choices of the top-squark mixing angle (9 t = 0.0, 0.10, 0.25, |). The dotted 
curve on the upper- left-hand panel corresponds to = m t ~ 2 = 50 GeV. Note the 
dips on the curves when m t pa 2 + m x . 
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Figure 6: The relative (%) one-loop supersymmetric electroweak-like correction to the 
top-quark pair-production cross section at the Tevatron as a function of the (higgsino- 
like) neutralino mass, for m t = 175 GeV, tan/5 = 1, = 75 GeV, m^ 2 = 250 GeV, 
and various choices of the top-squark mixing angle (9 t = 0.0, 0.10, 0.25, |). The dotted 
curve on the upper- left-hand panel corresponds to = m t ~ 2 = 75 GeV. Note the 
dips on the curves when m t pa 2 + m x . 
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Figure 7: The relative (%) one-loop supersymmetric QCD-like correction to the top- 
quark pair-production cross section at the Tevatron as a function of the universal 
squark mass, for m t = 175 GeV and the indicated choices of the gluino mass. The 
dashed curve represents the rrig = rriq case. 
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Figure 8: The relative (%) one-loop supersymmetric QCD-like correction to the top- 
quark pair-production cross section at the Tevatron as a function of the squark or 
gluino mass, for m t = 175 GeV and for the indicated choices of top-squark masses 
(raj , mf ). The dotted curves highlight the effect of non-degenerate top-squark 
masses. The dashed curve represents the fully degenerate 2 = rrig = rriq case. 



22 



